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ABSTRACT 


Theory  for  aeolian  tone  intensity  created  by  rigid 
circular  cylinders  has  been  found  to  agree  quite  well  with 
experimental  results.   This  investigation  is  to  experimentally 
explore  the  wake  of  a  cylinder  shedding  vertices  at  a  frequency 
that  is  equal  to  one  of  the  natural  vibrational  frequencies  of 
the  cylinder.   Investigation  was  carried  out  in  the  following 


areas: 


1,  lateral  spacial  correlation  of  vortex 
shedding, 

2,  distance  behind  cylinder  and  value  of 
RMS  velocity  at  assumed  position  where 
vortices  start  to  "roll  up", 

3,  observation  of  the  phase  betv;een  velocity 
signals  at  two  points  along  a  cylinder. 
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CHAPTER  I 
INTRODUCTION 

The  study  of  fluid  flow  past  circular  cylinders  has  been 
in  progress  since  the  latter  part  of  the  19th  century  when  Strouhal 
first  showed  that  the  vibrations  causing  the  sound  were  transverse 
to  the  direction  of  the  wind.   The  early  interest  in  the  phenomena 
was  mainly  concerned  with  the  vortex  shedding  frequency  as  a 
function  of  Reynolds  number. 

Until  recently,  very  little  work  has  been  done  concerning 
the  intensity  of  vortex  shedding  sounds  known  as  aeolian  tones, 
Phillips    and  Etkin,  Korbacher  and  Keefe    have  derived 

theoretical  expressions  for  the  radiated  sound  intensity  by  starting 

(3) 
with  Lighthills    basic  theory  of  aerodynamic  sound  as  had  already 

(M 

been  extended  to  rigid  bodies  by  Curie    . 

Phillips  used  the  product  of  £  £   to  account  for  the 

1  '  c 

fact  that  the  force  may  not  be  coherent  over  the  entire  length  of 
the  cylinder.   The  use  of  the  product  &  %     is  valid  if  the  body  is 
long  compared  to  the  coherent  length.   His  resulting  equation  for 
sound  intensity  is: 


— £~  S2  c7  '?   U3  l    *  sin28 
,  ,  2       *  c 


I    « 


(1) 


(-)  ~     indicated  time  averages 
p   ~     free  stream  density 
r,0   -  polar  coordinates  (9  being  the  upstream  direction). 

5  "  -r—  =  Strouhal  number 

L  o 

d  ■  cylinder  diameter 

N  =  vortex  shedding  frequency 

U  -  free  stream  velocity 

C.  c  lift  coefficient 

M  -  Mach  number 

6  -  length  of  cylinder  shedding  vortices 
£  =  spanwise  coherent  length  of  vortices 

I   -     acoustic  intensity 

Etkin,  Korbacher,  and  Keefe  derived  a  more  general  equation  that 
considers  any  length  of  coherent  force, 

In  order  to  allow  for  longer  lengths  of  coherent  force, 

the  centroid  of  the  area  under  the  correlation  coefficient  must 

2 
be  included.   The  effective  length,  that  is  substituted  for  I      is 

I      (£"Y)  rather  than  I   I    , 
c  c 

Y  "   centroid  of  the  area  under  the  cross  correlation 
coefficient  curve,  y   becomes  small  as  the  length  of  coherent 
force  becomes  small.   Therefore,  equation  (1)  is  apporached  for 


small   £    ,      The  more   general  equation  from  Etkin,   et   al»    shows: 


c 


I        »   _^»~   S2    C      H3  y3J,  (£_y)   Sin2e 

16  2       * 

(2) 

Experimental  investigation  of  radiated  sound  intensity 
from  circular  cylinders  when  shedding  vortex  wakes  has  been  made 
by  Gerrard    and  Keefe    ,   In  general >  these  experimental  re- 
sults confirm  the  theory  for  rigid  bodies  especially  in  the 
lower  Reynolds  number  region. 

The  flow  past  cylinders  that  are  non-rigid  bodies  has 
not  been  investigated  sufficiently  to  determine  what  happens  when 
the  vortex  shedding  frequency  and  a  natural  frequency  of  the 
vortex  producing  body  coincide.   This  thesis  is  a  report  of 
experimental  investigation  of  certain  properties  of  the  wake 
behind  a  cylinder  shedding  vortices  at  a  natural  vibrational 
frequency  of  the  cylinder  producing  the  vortices. 

The  theoretical  equation  mentioned  earlier  for  sound 

intensity  can  be  simplified  for  evaluation  to: 

I  -  C.2   (A-Y)A 
*>         c 

(3) 

which  shows  the  three  terms  that  may  change  for  the  situation 

of  frequency  coincidence  C«,  £  ,  and  Ye   Current  work  using  the 

*■   c 


same  experimental  set  up  is  proceeding  to  obtain  experimental  values 
for  sound  intensity,  lift  coefficient,  and  any  "locking  together" 
of  the  tv;o  frequencies,   Phillips    expresses  the  opinion  that  the 
correlation  length  may  increase  greatly  when  the  vortex  shedding 
frequency  and  the  vibration  frequency  coincide  due  to  a  coupling 
of  fluid  motions  at  different  axial  positions  caused  by  the  cylinder 
motion.   This  supposition  was  explored  experimentally  and  the  re- 
sults are  presented  here. 


CHAPTER  II 
PROCEDURE 


11,1,   Lateral  Spacial  Correlation 


The  aeolian  tone  test  rig,  which  is  pictured  in 
Appendix  A,  was  used  to  hold  a  cylinder  of  ,154  inches  diameter 
and  15  inches  in  length  under  very  high  tension  to  produce  the 
desired  natural  frequency.   It  was  desirable  to  have  one  of  the 
lover  harmonic  vibrational  frequency  coincide  with  the  vortex 
shedding  frequencies  in  order  to  eliminate  as  many  nodes  in  the 
cylinder's  vibration  as  possible.   The  cylinder  was  put  under  as 
high  tension  as  the  test  rig  was  physically  capable  of  holding. 
The  table  below  gives  the  natural  frequencies  of  vibration  and 
the  air  speed  and  Reynolds  number  required  to  cause  vortices  to 
shed  at  these  frequencies.   Only  the  wake  of  the  2nd,  and  3rd, 
harmonic  was  explored, 

D       Ud 


K  (cps)     U  (fps) 
2nd,  harmonic     1030       62,5       4,57  x  10 


e    v 

3 


3rd,  harmonic     1520        83        6.07  x  103 


The  cross-correlation  coefficient  was  computed  in  the 
following  manner: 


AIT  '  ,  4 


V?    V?     V^  ViT 


(4) 


R   (z)  B  the  cross-correlation  coefficient  as  a  function  of 


spacial  separation  of  the  hot  wire  anemometers. 

The  cross-correlation  coefficient  is  stated  here  to  be 
only  a  function  of  the  separation  and  not  a  function  of  where  the 
separation  is  located  along  the  cylinder.   The  investigation  to 
determine  if  this  coefficient  and,  in  turn,  the  correlation  length 
is  independent  of  location  along  the  cylinder  was  not  carried  out 
and  is  a  subject  for  future  investigation. 

The  values  of  R(z)  were  then  plotted  against  z.   See 
Figures  [4],  [5],  [6],  and  [7]«   A  mechanical  integrator  was  used 
to  determine  the  area  under  the  R(z)  curve.   This  value  is  defined 
as  one  half  the  correlation  length* 

i.e.    I     »   /    RAn(z)  dz 
c    J  Ab 

•00 

(5) 


A  plot  of  £  as  a  function  of  Reynolds  number  is  shown  on  Figure  [S], 
c 


In  this  range  of  Reynolds  numbers,  the  Strouhal  number  is  practi- 
cally constant  at  ,205,   The  cylinder  consisted  of  a  ,091"  diarreter 
steel  wire  over  which  was  placed  a  Teflon  tube  that  brought  the 
overall  diameter  to  ,154  inches. 

This  cylinder  was  placed  in  the  open  jet  test  section 
of  the  low  turbulence  wind  tunnel  in  the  M.I.T.  Acoustics  and 
Vibrations  Laboratory,   Figure  [1A] shows  the  dimensions  of  the 
wind  tunnel.   It  is  more  fully  described  in  a  report  by  Hanson    , 

For  the  correlation  measurements  one  of  the  hot  wire 
sensing  elements  was  placed  in  a  stationary  position  while  the 
other  sensing  element  was  held  on  a  traverse  mechanism  that  could 
be  moved  in  a  direction  parallel  to  the  cylinder.   The  two  hot 
wires  were  oriented  parallel  to  the  cylinder  and  positioned  at 
the  edge  of  the  wake  where  the  signal  is  periodic,  but  not 
turbulent.   This  position  varied  from  speed  to  speed,  but  in 
general  was  about  1/2  inch  back  and  1/4  inch  to  one  side  of  the 
center  of  the  cylinder.   See  Figure  [2],   The  signals  from  these 
two  Constant  Temperature  Anemometers  (CTA)  went  to  a  correlation 
unit  that  produced  a  sum  or  difference  output.   The  RMS  output  was 
recorded  on  a  Graphic  Level  Recorder,   Four  outputs  were  recorded 


VcaV  .  V^v  .  v? 


,2   A  ,2 


See  Figure  [3],   A,  B  =  voltage  signals  from  the  two  CTA's, 


To  fur tiier  check  the  correlation  at  coincidence  of 
frequencies  a  dual  beam  oscilloscope  was  used  to  observe  the 
phase  relation  of  the  signals  from  each  hot  wire. 

To  ensure  that  the  vortex  shedding  frequency  is 
coincident  with  a  natural  frequency,  a  sine  wave  signal  from  an 
oscillator,  set  on  the  desired  natural  frequency  by  an  electronic 
counter,  was  displayed  on  the  dual  beam  scope.   The  speed  of  the 
tunnel  was  then  adjusted  to  cause  the  hot  wire  signal  to  match 
that  frequency  of  the  oscillator. 

After  setting  the  wind  tunnel  speed,  the  signals  from 
both  CTA's  were  displayed  on  the  oscilloscope  for  observation  of 
the:  phase  between  the  two  signals. 


11,2,   Pes itionjgf ^  liaximum  RITS  Ve loc ity 

To  determine  the  position  where  it  is  assumed  that  the 
vortices  start  to  'roll  up"  one  hot  v;ire  anemometer  was  used  with  a 
two  directional  traverse  mechanism,   A  search  technique  was  used  in 
the  following  manner.   Starting  very  close  behind  the  cylinder ,  the 
wake  was  traversed  noting  the  position  at  which  the  maximum  RMS 
velocity  occurred.   The  value  of  RMS  voltage  and  the  DC  voltage 
output  of  the  anemometer  was  recorded  at  this  location  along  with 
the  position.   The  hot  wire  sensing  device  was  then  backed  away 
from  the  wire  a  little  further,   The  same  traverse  was  made 
recording  the  same  values  at  this  point.   This  was  repeated  until 
the  maximum  value  of  RMS  voltage  was  bracketed.   Then,  the  exact 
position  was  located  by  searching  within  the  bracketed  space,   A 
typical  plot  of  how  the  maximum  RMS  voltage  varies  as  a  function 
of  distance  behind  the  wire  is  presented  as  Figure  [9],   The  maxi- 
mum value  is  presumed  to  be  the  location  of  where  the  vortices 
start  "rolling  up", 

A  traverse  behind  the  cylinder  was  performed  at  three 
speeds  to  have  an  idea  of  the  properties  of  the  wake.   The  results 
of  RMS  velocity  and  mean  velocity  are  plotted  and  presented  as 
Figure  [10],  [11],  and  [12], 


CHAPTER  III 
DISCUSSION  OF  RESULTS 

III.].,  h^^S^.^S^^J^S^L^^Ji£^^2£Jfl. 

As  stated  in  the  introduction  the  equation  for  the  intensity 
of  aeolian  tones  is  proportional  to  lift  coefficient,  correlation  length,, 
and  centroid  of  the  area  under  the  cross-correlation  coefficient  curve 
as  thus: 


I  -  C^2  a-Y)*c 


(3) 


The  correlation  length  £  is  the  effective  length  of  a 
cohereiit  vortex.   As  the  Reynolds  number  increases,  the  two 
dimensional  qualities  of  the  vortex  shedding  start  to  disappear. 
In  other  words,  there  is  a  span-wise  variation  in  the  phase  of  the 
vortex  shedding*   The  average  coherent  length  of  shed  vortices  is 
used  in  calculating  the  aeolian  tone  intensity.   This  average  length 
is  obtained  by  a  procedure  using  two-point  correlation  of  velocity 
along  a  cylinder  as  described  in  Procedure, 

Correlation  of  velocity  variations  rather  than  force  were 
made  for  two  reasons.   First,  It  was  impractical  to  put  pressure 
taps  on  the  cylinder  due  to  the  small  diameter.   Secondly ,  because 
the  cylinder  has  an  induced  motion  the  pressure  tap  readings  would 
consist  of  a  component  due  to  the  motion  as  well  as  the  pressure 
caused  by  flow  fluctuations. 
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el  Baroudi    made  similar  cross-correlation  measurements  on 

a  rigid  cylinder.   His  data  ranged  in  Reynolds  number  from  about 

4  4 

10   to  4,5  x  10  ,   He  found  that  the  correlation  length  increased 

as  Reynolds  number  increased.   His  value  of  I     increasing  from 

3, 2d  at  Reynolds  number  104  to  6d  at  Reynolds  number  4,5  x  10  , 

•  .        Figure  [8],  shows  correlation  length  decreasing  with 

increasing  Reynolds  number  from  about  8d  at  Reynolds  number  of 

3  3 

6,5  x  10   to  5d  at  Reynolds  number  of  11  x  10  , 

No  apparent  change  in  the  trend  of  &     was  noted  at  points 
of  coincidence  of  the  vortex  shedding  frequency  and  a  natural  vibra- 
tional frequency  of  the  cylinder.   It  was  thought  that  with  the 

cylinder  vibrating  primarily  in  one  mode  that  the  coherent  length 

£  (9 ) 

would  increase  possibly  approaching  —   (n  -  mode  number) ,   Gerrard   , 

in  fact,  measured  the  correlation  coefficient  at  low  Reynolds  number 
and  found  it  to  be  of  the  following  shape: 


R(z) 
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In  the  measurements  obtained  during  this  investigation, 
the  value  of  the  correlation  coefficient  approached  zero  monotoni- 
cally  at  all  velocities  including  those  at  which  the  vortex  shedding 
frequencies  coincided  with  a  natural  frequency  of  the  cylinder. 

More  cross-correlation  measurements  should  be  made  vith 
particular  emphasis  on  having  the  vortex  shedding  frequency  exactly 
match  the  natural  frequency  and  also  more  measurements  at  frequen- 
cies near  the  natural  frequency. 

When  observing  the  signals  from  both  hot  wires  simultan- 
eously on  a  dual  beam  oscilloscope,  with  the  hot  wire  probes 
adjacent  to  one  another  the  signals  remained  inphase  as  would  be 
expected.   As  the  probes  were  moved  apart  parallel  to  the  cylinder, 
the  two  signals  did  not  remain  inphase  nor  show  a  constant  phase 
shift.   Instead  the  two  signals  would  move  in  and  out  of  phase  with 
no  apparent  pattern. 

To  check  for  phase  reversal  when  crossing  nodal  points 
the  hot  wires  were  placed  on  each  side  of  a  nodal  point  at  a  point 
midway  between  the  nodal  points  as  thus  for  the  3rd,  harmonic, 
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Again  the  signal  moved  continuously  in  and  out  of  phase 
as  previously  described. 

In  all  these  observations,  on  the  dual  beam  scope,  there 
didn't  appear  to  be  any  change,  of  observed  signal  whether  the 
vortex  shedding  frequency  was  in  coincidence  or  not  with  a 
natural  frequency. 
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III.  ,2.   Position  of  Maximum  RMS  Velocity 


The  position  of  the  maximum  RMS  velocity  is  assumed  to  be 
the  point  where  the  vortex  is  "rolling  up". 

The  maximum  RMS  voltage  referred  to  is  the  unfiltered  hot 
wire  anemometer  signal.   The  signal  at  the  position  of  maximum  RMS 
velocity  contained  random  frequencies  together  with  the  vortex 
shedding  frequency,   A  frequency  spectrum  was  taken  of  the  hot 
wire  anemometer  signal.   The  background  noise  was  at  least  15db 
lower  at  all  times  than  the  signal  at  the  vortex  shedding  frequency, 
Therefore^  the  addition  to  the.  broad  band  signal  was  less  than 
0,5db  or  ,01  volts.. 

As  shown  in  Figure  [9]  the  magnitude  of  the  fluctuating 
velocity  builds  up  linearly  as  the  hot  wire  is  moved  downstream 
to  a  maximum  then  falls  off  showing  the  viscous  diffusion  of  the 

vorticity. 

(9) 

Birkhoff    assumes  that  the  wake  acts  as  an  oscillating 

flat  plate  attached  to  the  cylinder  providing  alternating  lift  as 
the  wake  swings  from  side  to  side.  He  develops  a  satisfactory 
Strouhal  number  using  this  assumption.  If  it  is  assumed  that  the 
point  of  maximum  RMS  voltage  is  the  point  where  the  vortex  starts 
"rolling  up"  and  that  birkhoff 's  evaluation  of  the  wake  acting  as 
an  oscillator  providing  alternating  lift,  then  the  distance 
downstream  to  the  point  of  maximum  RMS  voltage  and  also  the 


14 


distance  off  to  one  side  of  centerline  where  the  maximum  RMS  voltage 
occurs  would  affect  the  lift  produced  by  the  wake. 

The  distance  downstream  to  the  point  of  maximum  RMS 
voltage  is  seen  from  Figure  [13]  to  decrease  linearly  in  the 
Reynolds  number  range  considered.   This  follows  the  same  trend  as 
reported  in  Goldstein     from  data  taken  by  Schiller  and  Linke, 

Again  we  see  no  deviation  from  this  trend  when  the  vortex 
shedding  frequency  and  a  nautral  frequency  coincide. 

Figures  [10,  [11],  and  [12]  show  the  mean  velocity  and 
RMS  velocity  of  the  wake  at  three  different  speeds  with  one  of  the 
situations  having  the  vortex  shedding  frequency  coincide  with  the 
3rd,  harmonic  of  the  natural  frequency.   The  results  show  that 
with  Reynolds  number  increasing  the  maximum  difference  in  mean 
velocity  increases  indicating  an  increase  in  vortex  strength. 

The  maximum  RMS  velocity  occurs  at  the  point  where  the 
mean  velocity  is  halfway  between  the  maximum  and  minimum  mean 
velocity.   This  is  expected  because  this  is  the  area  of  maximum 
shear  flow.   The  slope  of  the  mean  velocity  increases  in  the 
region  of  maximum  shear  flow  with  an  increase  in  Reynolds  number. 

In  conclusion,  there  did  not  appear  to  be  any  abnormal- 
ties  in  the  properties  of  the  wake  when  coincidence  of  frequencies 
occurred. 
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APPENDIX  B 
EQUIPMENT  USED  FOR  CORRELATION  LENGTH  MEASUREMENTS 

Thermo-Sys terns,  Inc, 

Model  1010A  Constant  Temperature  Anemometer 

Model  1015B  Correlator 

Model  1031-2  Monitor  and  Power  Supply  with 
(2)  Model  1035  Linearizer  modules, 

The  sensing  device  used  was  a  Tungsten  Hot  Wire  ,00015" 
diameter  and  ,05"  long. 

The  output  of  the  Constant  Temperature  Anemometer  was 
led  into  the  linearizer  modules  so  that  given  fluctuations  in 
velocity  would  give  the  same  voltage  f luctuations  for  any  mean 
velocity.   The  linearized  output  was  sent  to  the  correlator  unit 
that  summed  and  differenced  the  two  signals; also  the  separate 
signals  were  obtained  from  the  output  of  this  unit.   The  output 
from  the  correlator  was  sent  to  the  Brilel  and  Kjaer  Frequency 
Analyser  connected  to  the  Graphic  Level  Recorder  o:\   which  the 
signals  were  recorded  in  db  level.   See  Figure  [3],    The  output 
was  continuously  observed  on  an  oscilloscope.   This  was  converted 
to  RMS  voltage  and  used  In  the  cross-correlation  coefficient 
equation.   The  magnitude  of  fluctuating  velocity  was  not  computed 
as  the  conversion  factor  cancels  out  of  the  R(z)  equation. 


31 


The  free  stream  velocity  of  the  wind  tunnel  was  obtained 
by  placing  the  hot  wire  sensing  device  in  a  free  stream  portion 
of  the  open  test  section  arid  recording  the  output  voltage  from  the 
linearizer  to  use  on  the  previously  obtained  calibration  curve. 
From  this  velocity t    together  with  the  cylinder  diameter  and 
kinematic  velocity s  the  Reynolds  number  was  computed. 

The  method  used  to  check  for  coincidence  of  frequencies 
was  to  run  a  frequency  analysis  of  the  output  of  an  accelerometer 
mounted  on  a  support  of  the  wire.   When  the  two  frequencies 
coincided f  the  natural  frequency  would  be  reenforced  to  its 
highest  value.   A  more  accurate  method  of  matching  frequencies 
was  devised  described  on  page  (8) , 
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APPENDIX  C 
RECOMMENDATIONS 

1.  Take  correlation  coefficient  measurements  at  different  locations 
along  the  cylinder  to  determine  if  there  is  spacial  homogeniety. 

2.  Work  with  a  system  that  will  allow  a  greater  amplitude  of 
cylinder  vibration  to  determine  how  amplitude  of  vibration 
affects  the  correlation  length.   The  correlation  lengths  should 
increase  due  to  a  coupling  of  mechanical  and  hydrodynamic  forces, 
Tliis  could  be  accomplished  by  increasing  the  diameter  of  the 
cylinder  without  increasing  the  mass  very  much.   This  also 
would  allow  a  coincidence  of  frequency  at  the  1st,  harmonic, 

3.  Take  more  correlation  measurement  in  close  proximity  to  the 
point  of  coincidence  of  frequencies  to  determine  if  there  is 
any  deviation  from-  the  I     vs,  Reynolds  number  trend  observed 
in  other  regions  of  frequency, 

4.  Take  correlation  measurement  data  in  the  same  Reynolds  number 
range  that  el  Baroudi  did  to  see  what  changes  in  trends  or 
values  are  observed  for  a  non-rigid  cylinder. 
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